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Percolation theory is applied to the phase-transition dynamics of domain pattern formation in 
segregating binary Bose-Einstein condensates in quasi-two-dimensional systems. Our finite-size- 
scaling analysis shows that the percolation threshold of the initial domain pattern emerging from 
the dynamic instability is close to 0.5 for strongly repulsive condensates. The percolation probability 
is universally described with a scaling function when the probability is rescaled by the characteristic 
domain size in the dynamic scaling regime of the phase-ordering kinetics, independent of the inter¬ 
component interaction. It is revealed that an infinite domain wall sandwiched between percolating 
domains in the two condensates has an noninteger fractal dimension and keeps the scaling behavior 
during the dynamic scaling regime. This result seems to be in contrast to the argument that the 
dynamic scale invariance is violated in the presence of an infinite topological defect in numerical 
cosmology. 

PACS numbers: 67.85.Fg, 64.60.ah 47.37.+q, 98.80.Cq 


I. INTRODUCTION 

The phase-transition dynamics of spontaneous symme¬ 
try breaking (SSB) is an important problem discussed in 
different fields such as condensed matter physics, cosmol¬ 
ogy, and high-energy physics [ifi. It is believed that 
the SSB transition causes the nucleation of topological 
defects due to the Kibble-Zurek mechanism (KZM) in 
the early universe [El, 0 . Because it is difficult to exam¬ 
ine the KZM in our universe experimentally, the defect- 
nucleation problem has been attracting much attention 
in the systems of helium superfluids and atomic Bose- 
Einstein condensates (BECs), given the analogy in SSB 
between our universe and the quantum fluids II in. 

Linear and planar defects such as strings and domain 
walls form a complicated network in quenched phase- 
transition dynamics. The nucleated defect is destroyed if 
it is closed; e.g., a closed loop of a quantized vortex (cos¬ 
mic string) can shrink and finally collapse owing to the 
interactions with phonons and/or quasiparticles (gravi¬ 
tational waves) |13l-[lq . except for a specific case known 
as vortons [I|, [ 3 , [I^- An infinite defect, an open defect 
across the system, can survive. It is important to study 
infinite defects because they could survive for a long time 
after the SSB phase transition begins and thus can exert 
a dominant influence on the later dynamics. In addition, 
the solution to this problem can provide an answer to the 
interesting question of whether or not infinite defects are 
present in our current universe. 

In cosmology, it is suggested that percolation is im¬ 
portant for understanding the defect-survival problem 
[TsI, [l9|. Standard percolation (SP) theory [20| aims to 
describe the quite common problem of how a connected 
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element “percolates” (penetrates through the system). 
The critical scaling behavior appears when the occupa¬ 
tion rate p of the element considered approaches a thresh¬ 
old Pc over which there appears a percolating cluster. 
The critical behavior is characterized by a set of universal 
critical exponents, which describe e.g. the fractal behav¬ 
ior of the percolating cluster and the size distribution of 
non-percolating ones independent of microscopic details 
of the physical systems. Regarding a real scalar field as 
the order parameter due to SSB under consideration, a 
domain wall is formed along the boundary between two 
spatial domains with positive and negative values in the 
field. Therefore, an infinite domain wall can appear when 
positive and negative domains simultaneously percolate 
with Pc <0.5 since the simultaneous percolation does not 
occur with pc > 0.5. The study of infinite domain walls 
is fundamental to understanding infinite strings because 
an infinite string may be considered as the intersection 
of two infinite domain walls in the real and imaginary 
parts of a complex scalar field [20|. These arguments 
sound quite reasonable, and it is meaningful to tackle the 
defect-survival problem in such a direction by utilizing 
the quantum fluid systems. However, it seems difficult 
to treat the problem in later stages of phase-transition 
dynamics because of its complicated non-equilibrium de¬ 
velopment of order parameters. 

An empirical law based on the dy namic scaling hy¬ 
pothesis in phase-ordering kinetics sheds light on this 
problem. This law states that spatial patterns of order 
parameters in earlier stages of phase-transition dynam¬ 
ics are statistically similar to those in later stages if the 
patterns are scaled by the characteristic domain sizes. 
Since the theoretical analysis in SP theory is directly re¬ 
lated to the statistical property of spatial distributions 
of domains (clusters), the phase-ordering kinetics can be 
reinterpreted in terms of percolation problem. For the 
sake of convenience, let us refer to such a study as phase¬ 
ordering pereolation in this paper. 
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Within recent years Arenzon et al [22| studied the 
statistics of the areas enclosed by domain boundaries dur¬ 
ing the curvature-driven coarsening dynamics of a two- 
dimensional nonconserved scalar field by carrying out nu¬ 
merical simulations on the square-lattice ferromagnetic 
Ising model from a disordered initial state with p = 0.5. 
They showed that the time development of the distribu¬ 
tion of domain areas is described well by using a dynamic 
scaling function with a critical exponent of SP theory for 
Pc = 0.5. Additionally, it is surprising that results of 
SP theory can predict the final state of the Ising system; 
the existence probability of a percolating domain with a 
specified topology is related to the possibility for the sys¬ 
tem to develop into the metastable stripe state with open 
domain walls in the last stage of the relaxation dynamics 

Although the occupation rate p is not a conserved 
quantity in those studies lHlil , p is kept constant statis¬ 
tically during the time development due to the dynamic 
scaling law. Thereby, in the dynamic scaling regime 
where the system is described by the law, we expect that 
such problems of phase-ordering percolation would be ap¬ 
plicable in a similar manner to systems of a conserved 
scalar field, in which p is conserved exactly. Study of 
phase-ordering percolation for both conserved and non¬ 
conserved fields develops our understanding of SSB tran¬ 
sition dynamics to a higher level. However, the phase¬ 
ordering percolation for a conserved field is not well stud¬ 
ied. 

In this paper, we theoretically study phase-ordering 
percolation for a conserved field of density difference 
between two condensates in segregating atomic binary 
BECs. Only recently, the dynamic scaling law has been 
examined theoretically for domain pattern formations of 
the BEG systems in different situations [29143^ , and now 
we are better ready to pioneer the problem of phase¬ 
ordering percolation in the quantum fluids. By per¬ 
forming numerical experiments of segregating BECs, our 
finite-size-scaling (ESS) analysis revealed that the perco¬ 
lation threshold is Pc ~ 0-5 for the initial domain pattern 
caused by the dynamic instability in strongly segregated 
BECs, independently of the intercomponent interaction. 
The scaling behavior of SP theory of the initial pattern 
is sustained until the characteristic domain size grows in 
the order of the system size. We show that an infinite 
domain wall can exist as an interface between two per¬ 
colating domains in binary superfluids and found that 
the wall has a non-integer fractal dimension in a scal¬ 
ing regime during the phase-ordering development with 
p = 0.5. 

This paper is organized as follows. In Sec. II, we for¬ 
mulate this system and perform the ESS analysis based 
on SP theory for the initial domain patterns emerging 
from the dynamic instability. In Sec. Ill, the scaling be¬ 
haviors of the percolating domain and an infinite domain 
wall are described in a unified method independent of the 
intercomponent interaction. It is shown that the scaling 
behaviors are sustained during the phase-ordering devel¬ 


opment by rescaling the ESS analysis with the character¬ 
istic domain size. Section IV is devoted to the conclusion 
and discussion. 


II. INITIAL DOMAIN PATTERNS 
A. Formulation 

We consider a percolation problem of density domains 
of quasi-two-dimensional binary BECs in a square box. 
By introducing the order parameter {j = 

1,2) for the jth component at low temperatures, the sys¬ 
tem obeys the Gross-Pitaevskii (GP) Lagrangian 

£ = y* {ih'ipidt'ipi + — K — V) dxdy (1) 

with the kinetic energy density K = p 

and the potential energy density 

X. ^( 1 + 712 ). , x2 , ^( 1 - 712 ). x2 /o^ 

y= - ^-(ni+n2)^ +- - -(ni-712)7 (2) 

where m is the particle mass, and g and ^712 are the 
intracomponent and inter component coupling constants, 
respectively. 

The total energy j dxdy{K -\- V) and the norm Nj = 
f dxdyrij are conserved quantities during the time devel¬ 
opment. Eor 712 > 1, the second term on the right-hand 
side of Eq. m becomes negative with = rii — 772 7 ^ 0 , 
and binary BECs undergo a spontaneous breaking of 
the Z 2 symmetry with Ud = n or —77, leading to do¬ 
main formation [2514^. The density rij becomes al¬ 
most zero in the domains of the k{^ j)th. component. A 
straight domain wall is stabilized under the pressure bal¬ 
ance Pi = P 2 with the hydrostatic pressure Pj = gn‘jl2 
in the domain of the jth component. By neglecting 
the thickness and curvature of the walls, we may write 
77j ^ 77 = {NiPN2)/L‘^ at a place far from domain walls. 
Without loss of generality, we consider the first compo¬ 
nent as the percolation element discussed below; then, 
the occupation rate is defined as p = Ni/{Ni N 2 ). 

B. Characteristic scales 

Eigure (Ha) shows the characteristic domain patterns 
emerging during the early stages of the time develop¬ 
ment for different values of p obtained by numerically 
solving the GP equations derived from the Lagrangian 
O- We imposed the periodic boundary conditions in 
our numerical simulations for the sake of computational 
convenience. The initial state for the numerical compu¬ 
tation is the fully miscible state with nj{t = 0) = Nj/L‘^ 
with a system size L. In our simulation, the dynamic 
instability is triggered by initially introducing a random 
fluctuation with a small amplitude S. Eor a strong insta¬ 
bility, the initial domain pattern emerging in the early 
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FIG. 1: (Color online) (a) The initial domain patterns of the 
density difference Ud with I ^ lo for 712 = 2.0, L = 204.8^, 
and p = 0.475, 0.500, 0.525. Domain walls are represented by 
green regions with lower values of ~ 0. (b) The maximum 
domain of (a) is represented by a red (blue) region for the first 
(second) component. The maximum domain is defined as a 
domain with the largest area in each component in the nu¬ 
merical box without taking the periodic boundary condition 
into account. 


stages is characterized by the wave number kM of the 
Bogoliubov mode with the maximum imaginary part F m 
in the dispersion . 

For the case of (A^i—A^ 2 )/(^ 1 +^ 2 ) <C 1, we have kM ^ 
V 712 1 /C and Fm ^ (712 - l)/2r with ^ = %/^gmn 

and r = ?i/gn. Then, the instability is characterized by 
the length /q and time tq as follows: 


lo = 


To = 


\/7l2 - 1 

2'kt ^ n 

-In ^^77 

712 - 1 2(^2 


(3) 

(4) 


In fact, we found that the domain patterns become clear 
around t = tq with the mean interwall distance I ^ Iq 
[see examples in Fig. [Ha)]. Here, I is computed by nu¬ 
merically integrating the total length 


R = L^ll 


(5) 


of the domain-wall line, defined by = 0. Our numer¬ 
ical simulations were done with a square computational 
mesh. The domain-wall line of = 0 is defined as a col¬ 
lection of sides between neighboring meshes with > 0 
and < 0. A saddle point, an intersection of the lines 
in the mesh space, occasionally occurs when two domain 
walls are close to each other. We calculated the aver¬ 
age value fid of rid around the saddle point and then 
regard the point with > 0 (< 0) as a point occupied 
by the first (second) component connecting the diago¬ 
nal domains with rid > 0 (< 0)- For r ^ tq, l{t) obeys 
a characteristic power law, which is typical behavior in 
the kink formation dynamics observed in similar systems 
[ 29143 ^: however, this is not the main topic of our work. 


C. A relic wall between the maximum domains 

Here, we mention the main target—open domain walls 
that do not form closed loops. We are interested in long- 
range open walls, which are distributed across the sys¬ 
tem. A long-range open wall is realized as the inter¬ 
face between the two domains, both of which simultane¬ 
ously span the system in the same direction, i.e., from 
X = —00 to 00 (or from y = —oc to 00 ). Then, the 
open wall between the spanning domains is directed to 
the x{y) direction. We call such an open domain wall a 
directed open wall, and other open walls are referred to as 
nondirected open walls. A directed open wall is deformed 
into a straight wall owing to some dissipative mechanism, 
which smooths the domain-wall lines. To study the dy¬ 
namics of such an open wall is important since it could 
survive for long times and thus influences later devel¬ 
opment of the phase-transition dynamics. Closed walls 
shrink and subsequently vanish. A closed wall can be 
stabilized around a so-called coreless vortex, in which a 
component is trapped by the core of a quantized vortex 
in the other component. The core size of such a vortex 
is on the order of the thickness of a domain wall, and we 
found that not so many vortices appear in our system. 
Thus, the effect is not important to our analysis. 

A long-range open wall exists as an interface between 
the two domains, whose area is the largest in each com¬ 
ponent. In our dynamical system, evaluation of the sur¬ 
vival of open domain walls is not straightforward. A di¬ 
rected open wall can be transformed into a nondirected 
open wall after a few reconnections of the wall with other 
short walls. An example of such a nondirected wall is 
demonstrated in Fig. [T][b) for p = 0.500, where the wall 
sandwiched between the maximum domains runs from 
X = —Ll2 to y = I//2. Conversely, it is easy to create a 
directed wall from such a long nondirected wall. Thus, 
an open wall sandwiched between the maximum domains 
has the potential to survive as “a relic” of the transition. 
We call such a long-range wall the relic wall. We will 
evaluate the possibility of a relic wall in connection with 
the statistical properties of the maximum domains. SP 
theory is useful to the evaluation. 


D. Percolation probability 

To demonstrate the direct relation of our system with 
the percolation problem, we computed the ensemble av¬ 
erage of the area Sm{p^L) of the maximum domain. In 
SP theory [20], a percolation problem is analyzed by cal¬ 
culating the area of the maximum cluster divided by the 
system area (see also the Appendix). Then, our percola¬ 
tion problem is analyzed by the quantity 

P{p,L) = {SM{p,L))JLr (6) 

where (• • •)^ represents the ensemble average with I fixed. 
We took the ensemble average of all the domain pat¬ 
terns satisfying the condition \l{t) — lo\/lo ^ 0.005 during 
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FIG. 2: (Color online) (a) Plots of the probability P(p, L) for 
712 = 2.0 and L/^ = 25.6,51.2,102.4,204.8,409.6. (b) Plots 
of the parameter P(p, . The threshold pc is determined 

by averaging the positions of all intersection points between 
the fitting functions represented with bold solid lines. The 
inset is the finite-size-scaling plot of (a). 


each development by performing numerical simulations 
Nl times with Np > 6553.6 x ^/L. We found that the 
maximum domain spans the system probably if p ^ 0.5 
{p ^ 0.5) for the first (second) component [see Fig.[Tl[b)]. 
Actually, P{p^L) increases rapidly from p ^ 0.5. This 
upturn becomes sharper for a larger system size, and the 
probability asymptotically approaches P(p, L) ^ p for 
p —> 1. These behaviors are consistent with the finite- 
size effect in the SP systems [^ . 

To confirm the scalin g be havior of SP in our system, we 
perform a FSS analysis . It is convenient to utilize the 
property of the parameter P(p, with the critical 

exponents f3 = 5/36 and z/ = 4/3 [2^. This parameter 
takes a universal value independent of the system size L 
at p = Pc'. 

P{pc,L)L^/^ = const. (7) 

This property is used to estimate the percolation thresh¬ 
old [^. This estimation is sufficient for our main pur¬ 
pose, although the so-called Binder parameter is typically 
used to determine the critical point of the phase transi¬ 
tion . Figure [ 2 ] shows the p dependence of P(p, 
for several values of L. The plots with different values of 
L collapse onto a point at p ^ 0.5. 


E. Percolation threshold 

The percolation threshold is estimated by fitting a plot 
with a linear function in the least-squares method for 
0.5—< P < t).bp5pL with dpp = 0.005 x 409.6^/P [see 
in Fig.[2][b)]. We averaged the positions of all intersection 
points between the fitting lines with two different values 
of L and obtained Pc = 0.497(2) for 712 = 2.0. There 
are intrinsically large fluctuations around the threshold 
Pc ~ 0.5. Thus, to achieve the threshold value with satis¬ 
factory accuracy, we carried out many more simulations 
for — 5pL < p < 0.5 + We computed more than 


Np X 4 times for 0.485 < p < 0.515 and Np x 8 times for 
p = 0.5. 

According to the FSS hypothesis [ 2 ^, the quantity 
P{p,L) is scaled by a scaling function F{x) of x = 
Ap with Ap = p — Pc as 

P(p, P)P^/^ = P(P^/^Ap). ( 8 ) 

The scaling plot is successful when using the above value 
of Pc [see inset of Fig. [2][b)]. This result clearly shows 
that our system actually obeys SP theory. 

To gain a deeper insight, it is fruitful to investigate 
how the percolation behavior depends on the intercom¬ 
ponent interaction, which affects the correlation of the 
domain patterns. It is known that the percolation thresh¬ 
old is sensitive to the correlation between elements, i.e., 



where the authors investigated the probability of the si¬ 
multaneous percolation of the two components in the 
two-dimensional model of two-phase fluids and showed 
the percolation threshold can change depending on the 
correlation. 

In our system, namely, a two-phase superflow, the cor¬ 
relation is controllable via the intercomponent interac¬ 
tion: the coherence (or healing) length r\j 7 V 112 -1 
of the field Ud- However, we found that the percola¬ 
tion threshold remains Pc ~ 0.5 insensitively to 712 ; 
Pc = 0.498(4), 0.498(1), and 0.496(2) for 712 = 1.7, 2.5, 
and 3.0, respectively. The point is that segregated binary 
BECs are classified into strong segregation (712 ^ 1.7) 
and weak segregation (1 < 712 ^ 1.7), which are related 
to the competition between and the healing length of 
the total density ni + 712 [43 • Nevertheless, we did not 
find any qualitative change for 712 >1.7. 

It is believed that the percolation threshold of a 
symmetric random scalar field is Pc = 0.5 in a two- 
dimensional continuum according to the intuitive argu¬ 
ment in Ref. [39|. Our system is regarded as a contin¬ 
uum system since the numerical simulations were done 
with a square computational mesh whose size is small 
enough compared to any relevant length scale of this sys¬ 
tem. Supposing the field Ud plays the role of the random 
scalar field, our result seems to be consistent with the 
argument of Ref. . 


III. RESCALING ANALYSIS 
A. Rescaling plot 

The fact that the percolation property is not sensi¬ 
tive to the intercomponent interaction makes us expect 
that domain patterns with different 712 can be univer¬ 
sally described as having the same statistical property. 
The FSS analysis is applied after L is scaled by the char¬ 
acteristic domain size, that is the mean interwall distance 
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FIG. 3: (Color online) (a) The rescaling plot with 

the threshold pc- The probability P{p,L) with 
L/^ = 25.6,51.2,102.4,204.8,409.6 is rescaled for 

712 = 1.5,1.7, 2 . 0 ,2.5,3.0. (b) The L dependence of 

the relic probability Q (filled marks). Shapes of the filled 
marks corresponds to those of the legend symbol in (a) 
for reference of the values of 712 . A plot of the rescaled 
percolation probability P at p = 0.5 is displayed with the 
same marks as in (a). White crosses represent the dynamical 
plot of Q{L{t)) with 712 = 2.0 and L/^ = 204.8 for several 
values of L(t) = L/l{t). 

I = ^ 0 ( 712 )- Then, the probability in Eq. (j 6 ]) is rescaled 
by the effective system size L = L/l as 

P (p, i) = F (Zi/'^Ap) , (9) 

where we used the relation P(p, L) = P(p, L) and the 
rescaling function F. 

This conjecture is confirmed by the successful overlap¬ 
ping of the rescaling plots by using the average value 
of the thresholds of 712 = 1.7,2.0,2.5, and 3.0 and 
Pc ~ 0.497, in Fig. [3fa). Note that the numerical exper¬ 
iments become more difficult to perform with decreas¬ 
ing 712 because the computation time and the required 
system size are proportional to tq and /q, respectively. 
For this reason, we could not determine the percolation 
threshold for 712 = 1.5 with satisfactory accuracy. How¬ 
ever, our rescaling plot successfully includes the data for 
712 = 1.5. From this fact, we conclude that the perco¬ 
lation properties of this system are nearly independent 
of 712, although we cannot examine the limit 712 ^ 1 
because of the difficulty in numerical simulation with 
To ^ 00 and Iq 00 . 

B. Relic probability 

Now, let us return to the relic-wall problem. To evalu¬ 
ate the spatial distribution of the relic wall, we introduce 
the relic-wall length Rm and the probability that a spe¬ 
cific wall segment belongs to the relic wall: 

Q{L) = {Rm/R)i = I {Rm)i /LP ( 10 ) 

We consider the rate p = 0.5, with which the dynamical 
analysis is valid, as discussed later. Because this rate is 


very close to the percolation threshold, we assume that 
the two maximum domains have the scaling behavior de¬ 
scribed by SP theory. From Eqs. m and ([7]), we have 

Sm{p ^ VcL) <x. , ( 11 ) 

and the percolating domain has the fractal dimension 

Ds = 2- (12) 

Similarly, we assume the scaling behavior 

Rm{L)ocL^^ (13) 

with the fractal dimension Dp of the relic wall. 

In general, a fractal dimension D of an object in two- 
dimensional space takes values ranging from zero to two 
depending on the scale of view. For example, in our 
system, a closed domain wall with a finite size has D = 0 
if it is viewed on the scale of L ^ 00 . Similarly, if the 
domain wall forms a straight line across the system, we 
have D = 1. When a single domain wall is distributed 
uniformly throughout the system, its fractal dimension 
is D = 2. 

Because of the negligible deviation of the percolation 
threshold Pc from the critical rate (p = 0.5) of the phase¬ 
ordering kinetics, the fractal dimension Dp is restricted 
by the scaling behavior of the percolating domains. If 
domain walls are homogeneously distributed throughout 
the system, the relic probability Q{L) is considered to 
be the ratio of the area occupied by the relic wall to the 
system area L^. Because the relic wall is sandwiched 
between percolating domains, the area ratio must be 
smaller than the right-hand side of Eq. (©, which is the 
area ratio of the maximum domains to . Then, for 
large L, the following condition should be satisfied: 

1 < Dp < Ds. (14) 

The lowest limit represents the requirement Rm ^ T, 
coming from the fact that the length of a relic wall be¬ 
tween two percolating domains is not smaller than the 
linear size of the system. 

A rescaling analysis is convenient to produce a unified 
description of the scaling behavior of our system. From 
the fractal behavior of the maximum domain, we have 

P{p^Pc,L)^L-l^/\ (15) 

In the same way, the relic probability Q{L) = Q{L) is 
written as 


Q{L) - L-^ 

(16) 

with the exponent 


III 

to 

1 

b 

(17) 

Then, the condition (111)1 reduces to 


/3/z/ < h < 1. 

(18) 
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FIG. 4: (Color online) (a) A typical time evolution of L{t) 
and (b) the snapshots of the maximum domains in a nu¬ 
merical simulation for p = 0.5 at t/ro ~ 2,9, and 33 with 
L = 46.6,23.3, and 11.6, respectively, (b) The maximum 
domains are represented by a red (blue) region for the first 
(second) component in the same manner as Fig. [T](b). 


sion Dr of the relic wall in our system is larger than that 
(^ 1 . 6 ) of the longest disclination line in a quenched ne¬ 
matic liquid crystal [42|, where the cell of twisted liquid 
crystal is so thin that the system can be regarded as a 
two-dimensional system of a non-conserved field and the 
disclination as a one-dimensional domain wall . This 
discrepancy may be attributable to the two facts that the 
temperature in the cell gradually decreased during the 
relaxation dynamics and the system size was not large 
enough to neglect the fluctuation and to confirm clearly 
the scaling behavior. 


IV. CONCLUSION AND DISCUSSION 


If the relic wall does not exhibit any scaling behavior, Dr 
should be an integer, that is, = 1 from the restriction 
m with Dr = 1. 

Figure [3](b) shows the rescaling plots of Q{L) for the 
initial domain pattern with L = i^//o( 7 i 2 )- The condi¬ 
tion m is actually satisfied with h ~ 0.3. For com¬ 
parison, we plot P(p = 0.5, L), which is consistent with 
Eq. (p!5|) . Our scaling analysis is well applicable for 
L ^ 1. This is why P does not obey this law below 
L ^ 10. Our results show that a relic wall has scal¬ 
ing behavior with a noninteger fractal dimension with 
Dr ^ 1.7. 


C. Dynamic finite-size-scaling analysis 

The relic wall will maintain the scaling behavior for 
a long time during the phase-ordering development be¬ 
cause t he sy stem with p ^ 0.5 obeys the dynamic scaling 
law [U, [^ . The law states that the domain patterns at 
a later time are statistically similar to those at an early 
time if the patterns are scaled by the characteristic do¬ 
main size l{t). Then, we may replace L = L/Iq in the 
above analysis with the dynamical value L{t) = L/l(t). 

To test this argument, we computed the relic probabil¬ 
ity Q{L{t)) for /(t)//o = 0.8,1.0,1.4, 2.0, 2.8, 3.9, and 5.6 
by averaging many numerical experiments for 712 = 2.0 
and Lj^ = 204.8 [see Fig. [3][b) and Fig. 0]. We numer¬ 
ically confirmed that the relic wall maintains its scaling 
behavior with ^ 0.3 during the phase-ordering devel¬ 
opment until a later time with L(t) ^ 10. 

It is noteworthy that an infinite domain wall maintains 
its scaling behavior in our system; in contrast, the system 
of infinite cosmic strings is believed to not be scale invari¬ 
ant in numerical cosmology [iilli. There, the fractal 
dimension of the infinite cosmic string was estimated to 
be 2 approximately, connected with that of the random 
walk. In this sense, it is also interesting that the fractal 
dimension Dr is between those of the random walk and 
the self-avoiding walk, 2 and 4/3 [4l|, respectively. 

Finally, it should be mentioned that the fractal dimen- 


In conclusion, we found that the percolation threshold 
of segregating binary BECs is approximately 0.5. An infi¬ 
nite domain wall with a noninteger fractal dimension can 
survive in the dynamic scaling regime, where the char¬ 
acteristic domain size (/) is much larger than the system 
size (L). Experimentally, we will observe a relic domain 
wall in an atomic cloud whose size is larger than at least 
dozens of times the length ^ for 712 = 2.0. It is ex¬ 
pected that our scaling analysis is generally applicable 
to two-dimensional systems for both conserved and non- 
conserved fields. In the last stage with / < L, there 
appears a remarkable difference between the two fields. 
A large closed wall or a long-range open wall is left in 
the final state for conserved fields since p is fixed to 0.5, 
while the stripe states or the ground state with p = 0 or 
1 are realized for nonconserved field [ 2 ^. 

It is a nontrivial problem whether an infinite domain 
wall exhibits a noninteger-fractal behavior even when pc 
is not so close to 0.5. It must be instructive to investigate 
phase-ordering percolation in three dimensions, where 
the percolation threshold becomes smaller than that in 
two dimensions in general (2^. The application of our 
analysis to different phase-ordering systems would also 
be helpful to verify the exponent h and its relation with 
other critical exponents. 
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Appendix A: Critical behavior close to the 
percolation threshold 

Here, we make a brief introduction of the critical be¬ 
havior of the maximum cluster described by SP theory. 
The percolation probability is the probability that a given 
point belongs to an infinite cluster. In the limit of infinite 
system size L ^ oo, the probability is written as 

P{p)= (Al) 

where p and Sm are the occupation rate and the area of 
the maximum cluster of the element considered. 

According to the percolation theory, we have the fol¬ 
lowing critical behaviors around the critical point Pc^ 

P{p) - (p - Pc)^ (p > Pc), (A2) 

^p(p) ~ Ip-PcT'". (A3) 

Here, is the correlation length between two points in 
the same cluster. These equations reduces to 

P{p) ip > Pc)- (A4) 

The finite-size-scahng hypothesis states that the sta¬ 
tistical behavior in a system with a linear size L must 
be characterized by the ratio Lj^p. Correspondingly, we 
shall evaluate a quantity 

(A5) 

which exhibits the finite size effect through a scaling func¬ 
tion / around the critical point as 

P{p,L)^C/^"f{L/G)- (A6) 


From this formula, we suppose a scaling law 

P(p, = F{L^/''/^p) (A7) 

with Ap = p — Pc- 

The system size L is only a relevant length scale at 
p = Pc since diverges as (x \p — Pc\~^ with z/ > 0. 
Thus, P{pc^L) should be a function of not ^p but L, and 
we have 

/(A/Cp) ~ (A/^p)“'^A (p ^ (A8) 

By substituting Eq. (IMll with the above formula into 
Eq. (IA5p . we have 

Sm (p ^ Pc, A) - , Ds = 2- /3/p (A9) 

Here, this equation shows that the maximum cluster at 
the critical point has a fractal dimension Ds- 

To derive the formula of Eq. ©, we introduce 
proportional constants Cp and Cf as P{p,L) ^ 

Cpip'GGlip) and /(i/^p) ^ for 

p ^ Pc- Erom these relations one obtains 

P{pc^L)L^^^ = const. (AlO) 

This shows that the quantity P(p, should ap¬ 

proach a universal value at p = Pc independent of L. 
This argument is consistent with the scaling law dSB 
with Ap = 0. 
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